Number fluctuations of sparse quasiparticles in a superconductor 
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We have directly measured quasiparticle number fluctuations in a thin fllm superconducting Al 
resonator in thermal equilibrium. The spectrum of these fluctuations provides a measure of both the 
density and the lifetime of the quasiparticles. We observe that the quasiparticle density decreases 
exponentially with decreasing temperature, as theoretically predicted, but saturates below 160 mK 
to 25-55 /fim^. We show that this saturation is consistent with the measured saturation in the 
quasiparticle lifetime, which also explains similar observations in qubit decoherence times. 



In a superconductor the density of unpaired electrons 
(quasiparticles) should vanish when approaching zero 
temperature [T]. This crucial property promises long de- 
coherence times for superconducting qubits '2] and long 
relaxation times for highly sensitive radiation detectors 
[3] . However, relaxation times for resonators [H [S] and 
qubit decoherence times were shown to saturate at 
low temperature. Recent modeling [51 [3] suggests that 
non-equilibrium quasiparticles are the main candidate for 
this saturation, which was tested qualitatively by inject- 
ing quasiparticles into a qubit [10 . A direct measurement 
of the number of quasiparticles and the energy decay rate 
in equilibrium at low temperatures would provide new in- 
sight in superconductivity at low temperatures, crucially 
needed in the aforementioned fields. 

At finite temperature, it follows from thermodynamics 
that the density of quasiparticles fluctuates around an 
average value that increases exponentially with temper- 
ature |11| . Here we report a measurement of the spec- 
trum of these fluctuations in a single aluminium super- 
conducting fllm {Tc — 1.1 K) in equilibrium, for temper- 
atures from 300 mK to 100 mK. The number fluctuations 
show up as fluctuations in the complex conductivity of 
the fllm, probed with a microwave resonator. The spec- 
trum of these fluctuations provides a direct measure of 
the number of quasiparticles in the superconductor. We 
observe that the quasiparticle density decreases exponen- 
tially with decreasing temperature until it saturates at 
25-55 /im~'^ below 160 mK. We prove that the measured 
saturation of the quasiparticle lifetime to 2.2 ms below 
160 mK is consistent with the saturation in quasiparticle 
density. Besides the fundamental signiflcance, our exper- 
iment shows that it is possible to reach the fundamen- 
tal generation-recombination noise limit in microwave ki- 
netic inductance detectors based upon Al resonators. 

In a superconductor in thermal equilibrium, the den- 
sity of quasiparticles per unit volume is given by 



= 2Noy/2TrkBTAexp{-A/kBT), 



(1) 



valid at /cbT < A, with A^o the single spin density of 
states at the Fermi level (1.72 x 10^° ^m^^eV"^ for Al), 
ks Boltzmann's constant, T the temperature and A the 



energy gap of the superconductor. Two quasiparticles 
with opposite spins and momenta can be generated from 
a Cooper pair by a phonon with an energy larger than 
the energy gap. When two quasiparticles recombine into 
a Cooper pair, a phonon is emitted. These processes, 
schematically depicted in Fig. [T^, are random processes 
in equilibrium. Assuming a thermal distribution of quasi- 
particles and phonons at low temperature, the average 
quasiparticle recombination time is given by |12) 
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where Tc is the critical temperature of the superconduc- 
tor and To a material dependent, characteristic electron- 
phonon interaction time. Eqs. [l] and [2] predict a very 
low quasiparticle density and consequently a very long 
quasiparticle lifetime at temperatures T < Tc/10. 

The process of random generation and recombination 
of charge carriers is a well-studied phenomenon in solid 
state physics, in particular in semiconductors, but has 
hardly been studied in superconductors. In one ear- 
lier experiment, this generation-recombination noise was 
identified in the current fluctuations through a tunnel 
barrier connected to a small Al quasiparticle box [13j . 
although only down to an intermediate 210 mK. The 
general theory of quasiparticle number fluctuations in a 
superconductor [111 113] connects the frequency depen- 
dence of the fluctuations to the microscopic dynamics of 
quasiparticle generation and recombination. The dom- 
inant timescale of these processes is the recombination 
time of a quasiparticle, (r^ is about 1 ms in an Al film 
[i|), because the phonon pair breaking time and the 
phonon escape time are both much shorter, about 10"^*^ 
s based on a film thickness of 40 nm [12j [14] . In thermal 
equilibrium, the generation and recombination rates are 
equal and the variance of the random number fiuctua- 
tions cr^ =< SN'^p >= Nqp = UqpV, with V the volume 
of the system. The power spectral density of these fluc- 
tuations shows a Lorentzian spectrum, given by 
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FIG. 1. (color online) (a) Schematic of the process of gen- 
eration (left) and recombination (right) of quasiparticles in a 
superconductor. The single arrow symbolises a quasiparticle, 
the double arrow a Cooper pair and the wavy arrow a phonon. 
(b) A microscope picture of the microwave resonator, two 
times reduced in length for visibility. The open ends set the 
half wavelength resonance condition. The coplanar waveguide 
(CPW) through line is used for excitation and readout (from 
contact 1 to 2) of the resonator with a microwave signal, (c) 
The real and imaginary parts of the complex transmission S21 
as a function of frequency around the resonator resonant fre- 
quency. The gray arrow indicates the direction of increasing 
frequency. The 5*21 traces out a circle in the complex plane, 
off-set from the origin, around the resonant frequency of the 
resonator. We define a resonator amplitude A and phase 9 
with respect to the resonance circle centre as indicated. In 
red, the measured fluctuations in S21 are shown in case the 
readout frequency is fixed at the resonant frequency. The fluc- 
tuations far off-resonance are shown in blue, which are used 
as a calibration for system noise contributions. 

with u! the angular frequency. Eqs. [l]and|2]show that the 
product NqpTr is constant over temperature, whereas the 
total integrated power spectral density increases expo- 
nentially with temperature. This is because the band- 
width of the fluctuations increases exponentially with 
temperature as well, as it scales with . We emphasise 
that this property is unique for quasiparticle generation- 
recombination noise in a superconductor. 

We measure the quasiparticle number fluctuations us- 
ing a high-quality microwave resonator. The high fre- 
quency response of the superconductor is controlled by 
the quasiparticle density through the complex conduc- 
tivity (Ti — i(T2. The real part, cri, is resistive and de- 
notes the conductivity by quasiparticles. The imaginary 
part, CT2, is inductive and due to the superconducting 
condensate, the Cooper pairs |15j . Quasiparticle number 
fluctuations will show up as fluctuations in the complex 
conductivity. To measure the complex conductivity, a 
40 nm thick Al film was patterned into microwave res- 
onators. The film was sputtered deposited onto a C-plane 
sapphire substrate. The critical temperature is 1.11 K, 
from which the energy gap A = l.TGfcsTc = 168 fj,eV. 
The low temperature resistivity p — 0.8 /xf2cm and the 



residual resistance ratio RRR = 5.2. The film was pat- 
terned by wet etching into distributed, half wavelength, 
coplanar waveguide resonators, with a defined central line 
width of 3.0 ^m and gaps of 2.0 ^m wide (Fig. [T]d). 
The resonator under consideration shows its lowest order 
resonance at 6.61924 GHz and has a central strip vol- 
ume of 1.0 • 10^ iim^. The low temperature (100 mK) 
resonance curve shows a coupling limited quality factor 
of 3.87x10^. The samples are cooled in a pulse tube 
pre-cooled adiabatic demagnetization refrigerator. The 
cold stage is surrounded by a superconducting magnetic 
shield inside a cryoperm shield. Special care has been 
taken to make the setup light tight, such that no excess 
quasiparticles are created by stray light. The sample is 
mounted inside a light tight holder which itself is placed 
inside another light tight box (also at base temperature) , 
to prevent radiation leaking in via the coax cable con- 
nectors. Radiation absorber, consisting of a mixture of 
carbon black, epoxy and SiC grains is placed inside both 
the sample holder and the outer box. The outer box is 
equipped with special coax cable filters that attenuate 
all frequencies above 10 GHz exponentially (see Ref. IH] 
for details). The system is proven to be light tight by 
measuring the quasiparticle lifetime as a function of the 
4 K-stage temperature while keeping the sample at base 
temperature. Within the measurement accuracy (15%), 
there was no change in the lifetime, indicating that the 
stray-light power at the chip is negligible. For the life- 
time measurements the sample is illuminated with a short 
pulse of light from a GaAsP LED (1.9 eV), as described 
in Ref. jH The fibre coupling of the LED to the sample 
is done via a 0.15 mm diameter waveguide that acts as 
a 1 THz high pass filter, to prevent any pair breaking 
radiation from the 4 K thermal environment where the 
LED is mounted to reach the sample holder [17] . 

The complex transmission of the microwave circuit is 
measured with a quadrature mixer and traces out a cir- 
cle in the complex plane. The microwave signal is am- 
plified at 4 K with a high electron mobility transistor 
(HEMT) amplifier and with a room temperature ampli- 
fier, before it is mixed with a copy of the original sig- 
nal [31 [T5] . We define a resonator amplitude and phase 
with respect to the resonance circle, as depicted in Fig. 
[1];. The resonator amplitude predominantly responds to 
changes in ai [H] (therefore also called dissipation di- 
rection). The responsivity of the resonator amplitude to 
quasiparticles dA/dNqp = —2aQK/V , with Q the quality 
factor of the resonator, a the fraction of kinetic induc- 
tance over the total inductance and k — ^^^^l^ which 
depends only weakly on temperature. The amplitude re- 
sponsivity was determined experimentally as described 
in Ref. HSl dA/dNqp is measured to be almost tempera- 
ture independent. For similar resonators it is known that 
the sensitivity in phase is limited by two-level fiuctuators 
[201 121] and that the sensitivity in amplitude is up to a 
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factor 10 better, limited by the HEMT amplifier [15] . 

The power spectral density due to quasiparticle num- 
ber fluctuations in the resonator amplitude is given by 
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where Tres is the resonator ringtime given by r^es — 
In this experiment 3> Tres ~ 2 fj,s, meaning that the 
roll-off in the noise spectrum will be determined solely 
by Tr if Sa is dominated by quasiparticle number fluc- 
tuations. Using Eqs. [T]|4j with tq = 438 ns [H] and a 
measured dA/dNqp = 5.0 x 10"'', we expect Sa = —99.3 
dBc/Hz. This is a high value compared to other super- 
conductors like Ta and Nb, due to the large tq in Al [T^ . 

We have measured the fluctuations in the resonator 
amplitude in equilibrium at the resonant frequency us- 
ing a microwave power of -77 dBm. The power spectral 
density, corrected for system noise, is shown in Fig. [2^ for 
various temperatures, which is the central result of this 
paper. In Fig. [2|d the system noise spectrum is shown, 
which is subtracted from the spectrum measured on res- 
onance, to get the corrected spectra in Fig. [2^. Parts of 
the time domain signal where large energy impacts are 
observed, are removed from the analysis as shown in Fig. 
[2];, because they distort the dynamic equilibrium. These 
impacts do not appear in the system noise signal and de- 
cay within a time r^, which shows they are due to events 
that create quasiparticles, for example due to cosmic ray 
hits [22] or local radioactivity. 

In Fig. [2ji we observe that the measured power spectral 
density of the fluctuations has a constant level for all 
temperatures. The roll-off in the spectra can be described 
with a single timescale that decreases with temperature. 
From these two properties we conclude that we directly 
observe quasiparticle number fluctuations. 

The measured recombination time, r^, is extracted 
from the resonator amplitude noise spectra and shown 
as the black squares in Fig. [3] as a function of temper- 
ature. At temperatures from 180-300 mK, we find the 
expected exponential temperature dependence. Eq. [2] 
is used to fit for the characteristic electron-phonon in- 
teraction time To and we find a value of 458±10 ns, in 
reasonable agreement with other studies [TU [23j . Due to 
the phonon trapping effect, which we cannot estimate 
accurately, the measured tq may differ from the pure 
electron-phonon time [TH [53]. At temperatures < 150 
mK we measure a temperature independent quasiparticle 
recombination time of 2.2 ms, which is among the longest 
reported for a thin superconducting film. Alternatively, 
the recombination time is measured by monitoring the 
restoration of equilibrium after a short pulse of optical 
photons [T , which is shown by the filled circles in Fig. [s] 
The lifetimes obtained from the noise spectra are equal to 
the lifetimes from the pulse measurement up to 220 mK, 
indicating that both measurements really probe quasi- 
particle recombination in equilibrium. The lifetime from 
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FIG. 2. (color online) (a) Corrected power spectral density 
of the resonator amplitude fluctuations as a function of fre- 
quency, calibrated for system noise contributions as depicted 
in (b). The spectra are plotted for six different temperatures. 
Note that the spectral density is temperature independent 
up to the roll-off frequency. The Lorentzian fits to the spec- 
tra, plotted as dashed lines, show that the spectra can be 
described with a single timescale. (b) The noise spectrum on 
resonance is corrected for system noise contributions by sub- 
tracting the system noise contributions from the raw data. 
The system noise is obtained by taking a calibration mea- 
surement at a frequency far from the resonance frequency as 
indicated in Fig. [ij;. (c) Part of the time trace of the resonator 
amplitude. High energy impacts are observed every 20-30 s, 
creating milUons of quasiparticles in the resonator volume. 
Smaller impacts, as shown in the inset, happen every 5-10 s. 
The time domain trace is filtered with a moving average filter 
with a time constant Tr/2, to distinguish small impacts from 
the noise. Parts of the time domain trace with impacts larger 
than 5 times the standard deviation are rejected, correspond- 
ing to energies of 0.43 eV or 2500 quasiparticles. 



the noise spectra agrees well with theory (full line) up 
to 300 mK. The lifetime obtained with the pulse method 
shows a deviation from theory, which we always observe 
in Al on sapphire resonators. The deviation between the 
two sets of experimental data calls thus, most likely, for 
a future analysis of the physical processes in the pulse 
method. 

We combine the level of the power spectral density 
(Fig. [2]) and the quasiparticle lifetime obtained from the 
roll-off in this spectral density (Fig. |3]) to obtain the 
number of quasiparticles, Nqp, by using Eqs. [3]and|4) 
together with the measured values of dA/dNqp. The re- 
sult is plotted in Fig. |4] with the black squares. As a 
cross check, we convert the quasiparticle lifetime mea- 
sured from the noise roll-off directly into quasiparticle 
number by using Eq. [2] which is shown in Fig. [4] with 
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FIG. 3. (color online) The quasiparticle lifetime is obtained 
from the roll-off frequency of the resonator amplitude noise 
spectrum (Fig. |2| and plotted as a function of temperature. 
The solid line is the lifetime, calculated from theory. Addi- 
tionally the lifetime is determined by measuring the decay of 
the amplitude while restoring the equilibrium after a short 
light pulse. The exponential decay of the excitation is fit and 
the obtained decay time is plotted as a function of tempera- 
ture. The measurements show consistently that the lifetime 
saturates to about 2.2 ms below 160 mK. The error bars rep- 
resent statistical uncertainties obtained from the fitting pro- 
cedures. 



red triangles. We assume that the relevant volume is the 
central strip volume of the resonator. The quasiparti- 
cle number, obtained via these two methods consistently 
shows a saturation, giving a low temperature quasiparti- 
cle density of 25-55 ^m~^. We conclude that the quasi- 
particle lifetime saturation is due to a saturation in the 
quasiparticle density, consistent with the conjecture of 
Martinis et al. [8^. 

The question remains what the source of the non- 
thermal (in view of Eq. [I]) quasiparticle density below 
160 mK is. In literature, it is usually referred to as non- 
equilibrium quasiparticle density, which is inferred from 
saturating tunnel rates O [53] and attributed to electro- 
magnetic noise [5SJ [SS] and radiation [57] or, if the first 
two are eliminated, to cosmic rays, local radioactivity, 
slow heat release or stray light jS]- In our experiment, 
excess quasiparticles due to cosmic ray hits are excluded 
as shown in Fig. [2];. Local radioactivity is also excluded 
as far as it generates similar distinguishable events. The 
sample box and coaxial cables are thoroughly shielded 
from stray light, validated by the observation that is 
independent of the 4 K-stage temperature, as explained 
before. Additionally, if there would still be stray-light, it 
would cause an additional photon shot noise [3] contribu- 
tion, which decreases with increasing temperature, due to 
the decreasing quasiparticle lifetime. The observed tem- 
perature independent noise level therefore excludes stray 
light. We cannot completely exclude quasiparticle gen- 
eration by the microwave signal. The microwave power 
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FIG. 4. (color online) The number of quasiparticles as 
determined from the noise level as a function of tempera- 
ture. The inset shows that the noise level (Sa) divided by 
the responsivity is almost temperature independent. The er- 
ror bars indicate combined statistical uncertainties. Alter- 
natively, the number of quasiparticles in the superconductor 
is determined from the measured quasiparticle lifetimes. We 
observe that the number of quasiparticles saturates at around 
30,000-50,000 at low temperature, where an exponential de- 
crease is expected, as shown by the solid line. 



range over which we can measure the quasiparticle fluctu- 
ations is only about 4 dB due to limitations in resonator 
power handling and amplifier noise temperature. Over 
this power range the quasiparticle density is power inde- 
pendent. Intriguingly, the remnant quasiparticle density 
is so low that a rather unexplored regime of large spatial 
separation between quasiparticles is entered. 

For superconducting radiation detectors, the noise 
equivalent power (NEP) is a common way to express the 
sensitivity. From the noise level and lifetime measure- 
ments we determine [3] an electrical NEP of 3.3x10"^^ 
WHz~^/^ due to generation-recombination noise, not 
taken into account the system noise, which will increase 
the NEP by about a factor of two. For Al resonators with 
this geometry and the observed remnant quasiparticle 
density, this is the fundamental limit to the sensitivity. 
A significant improvement in the NEP is only possible 
if one could reduce the remnant quasiparticle density or 
the resonator volume. 
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